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Abstract 

Diffractive dissociation of virtual photons, 7*p — > Xp, has been studied in ep in- 
teractions with the ZEUS detector at HERA. The data cover photon virtualities 
0.17 < Q 2 < 0.70 GeV 2 and 3 < Q 2 < 80 GeV 2 with 3 < M x < 38 GeV, where M x 
is the mass of the hadronic final state. Diffractive events were selected by two meth- 
ods: the first required the detection of the scattered proton in the ZEUS leading 
proton spectrometer (LPS); the second was based on the distribution of Mx- The 
integrated luminosities of the low- and high-Q 2 samples used in the LPS-based anal- 
ysis are ~ 0.9 pb _1 and ~ 3.3 pb _1 , respectively. The sample used for the Mx-based 
analysis corresponds to an integrated luminosity of ~ 6.2 pb _1 . The dependence of 
the diffractive cross section on W, the virtual photon-proton centre-of-mass energy, 
and on Q 2 is studied. In the low-Q 2 range, the energy dependence is compatible with 
Regge theory and is used to determine the intercept of the Pomeron trajectory. The 
W dependence of the diffractive cross section exhibits no significant change from the 
low-Q 2 to the high-Q 2 region. In the low-Q 2 range, little Q 2 dependence is found, a 
significantly different behaviour from the rapidly falling cross section measured for 
Q 2 > 3 GeV 2 . The ratio of the diffractive to the virtual photon-proton total cross 
section is studied as a function of W and Q 2 . Comparisons are made with a model 
based on perturbative QCD. 
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1 Introduction 



The properties of high-energy hadron-hadron cross sections, notably the energy depen- 
dence of the total and elastic cross sections, are described successfully by Regge phe- 
nomenology in terms of the exchange of the Pomeron trajectory, ap{t) = op(0) + a'jpt, 
where t is the squared four- momentum carried by the exchange [jj]. The intercept and 
slope of the trajectory were found to be op(0) = 1.08 and a'p = 0.25 GeV" 2 , respectively, 
by Donnachie and Landshoff [fj| using the energy dependence of the hadron-hadron total 
and elastic cross sections. Such a Pomeron trajectory is referred to as "the soft Pomeron". 
At high energies, hadron-hadron total cross sections, including the jp total cross section, 
can be expressed in terms of this trajectory as 

where W is the virtual photon-proton centre-of-mass energy. 

Measurements of the diffractive dissociation of photons have shown that, for quasi-real 
photons (Q 2 ~ 0, photoproduction, where Q 2 is the exchanged photon virtuality), the 
value of qjjp(O) is compatible with the expectations based on soft-Pomeron exchange ||, 
U). The study of diffractive processes in ep collisions at large virtualities has opened 
up the possibility of investigating the Pomeron in a regime where perturbative QCD 
(pQCD) is applicable ||. In this regime, the exchange of the Pomeron trajectory may 
be described, at lowest order, as two-gluon exchange in the t channel so that the cross 
section is proportional to the square of the gluon density in the proton. Since the gluon 
distribution rises steeply at small Bjorken x (or, equivalently, for large values of W), a 
possible signature of the transition from the soft non-perturbative regime to the hard 
perturbative regime is a change to a M^-dependence of the cross section steeper than that 
from the exchange of a soft-Pomeron trajectory. The value of the Pomeron intercept, 
Ofp(0), measured in the deep inelastic scattering (DIS) regime (Q 2 ^ a few GeV 2 ) is larger 
than that of the soft Pomeron ||[7|], which suggests that pQCD effects have become 
important. 

In analogy with the usual DIS formalism for the proton structure function, F 2 , one can 
introduce a diffractive structure function, F 2 D . Studies of photon diffractive-dissociation 
have shown that, for Q 2 <; 1 GeV 2 , F 2 D has only a weak, logarithmic, dependence on 
Q 2 H0. However, conservation of the electromagnetic current requires that both F 2 
and F 2 must behave like Q 2 as Q 2 — > 0. 

In this paper, the inclusive diffractive dissociation of virtual photons, 7*p — > Xp, is 
investigated by studying the reaction ep — > eXp at HERA both in the perturbative region 
(Q 2 3> 1 GeV 2 ) and in the transition region between the non-perturbative (Q 2 ~ 0) and 
perturbative regions. The measurements are presented as a function of W and Q 2 . The 
Pomeron intercept is determined through the measurement of the energy dependence 
of the diffractive cross section in the transition region, which has not previously been 
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explored in diffraction. The W and Q 2 behaviour of the diffractive cross section and of 
the virtual photon-proton total cross section, o"7 ot p , are compared by studying their ratio 
as a function of W and Q 2 . 

Diffractive events were selected by two methods. The first required the detection of the 
scattered proton in the ZEUS leading proton spectrometer (LPS) and is referred to as 
the "LPS method" . Although statistically limited because of the small acceptance of the 
LPS, this method permits the selection of events with negligible background from the 
double-dissociative reaction, ep — > eXN, where the proton also diffractively dissociates 
into a state N of mass Mx that escapes undetected in the beam pipe. The LPS method 
also gives access to higher values of Mx, the mass of the hadronic final-state system, X, 
and allows the measurement of the squared four-momentum transfer at the proton vertex, 
t. The second method, henceforth referred to as the U M\ method" [0], is based on the 
characteristics of the distribution of Mx- The sample selected with the M x method 
contains a background contribution from the double-dissociative events. 

This paper presents results in the region 0.17 < Q 2 < 0.70 GeV 2 , obtained with both 
methods, and in the region 3 < Q 2 < 80 GeV 2 using only the LPS. The measurements 
cover the region 3 < Mx < 38 GeV. Results in the DIS region obtained using the M\ 
method have been previously reported 0. 



2 Kinematic variables and cross sections 

Inclusive diffractive dissociation of virtual photons in positron-proton collisions, ep — > 
eXp, can be described by the kinematic variables Q 2 , W, Mx, and t. The differential 
cross section for j*P Xp is related to the cross section for the reaction ep — > eXp by || 

d^Z(Q 2 ,W,M x ,t) = w) d 2 ^(Q 2 ,W,M x ,t) 
d\nQ 2 d\nWdM x dt W ' ' dM x dt ' 1 ' 

where 

r = -[i + (i-y) 2 ] 

is the virtual photon flux, a is the fine-structure constant, y ~ (W 2 + Q 2 )/s is the fraction 
of the positron energy transferred to the proton in its rest frame, and s is the square of 
the positron-proton centre-of-mass energy. 

In analogy with the formalism of inclusive deep inelastic ep scattering, the diffractive cross 
section for the reaction ep —>■ eXp can also be expressed in terms of diffractive structure 
functions 0: 

d*aZ _ Ana 2 f y 2 
d{3dQ 2 dx p dt (3Q 4 \ U 2(1 + R D W((3,Q 2 ,x F ,t)) 



,D(4) 



((3,Q 2 ,x p ,t), (2) 
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where the diffractive structure function and the ratio of the cross sections for lon- 

gitudinal and transverse photons, R D ^\ have been introduced. 

The variables x p and (3 are related to Q 2 , W 2 , M\ and t by 

Q 2 + M 2 x -t 



Q 2 + W 2 - M 2 

Q 2 

Q 2 + M 2 X -V 



where M p is the proton mass. The variables x p and (3 can be interpreted, assuming the 
t-channel exchange of a Pomeron with partonic structure, as the fraction of the proton 
momentum carried by the Pomeron and the fraction of the Pomeron momentum carried 
by the struck parton, respectively. 

Equations flip and (g) can be combined to give 



dM x dt Q 2 + W 2 Q 2 + M 2 Q 



where \t\ <Q 2 + M| has been assumed and M p and R D ^ have been neglected [10]. An 



analogous expression holds for the three-fold differential diffractive structure function, 
F^ 3 \ obtained by integrating over t. Equation (||D is the diffractive analogue of 

the expression o"7 ot p = (4:7r 2 a/Q 2 )F 2 which holds for inclusive 7*p scattering at high W. 



3 Experimental set-up 

The measurements were performed at the HERA ep collider at DESY between 1995 and 
1997 using the ZEUS detector. At that time, HERA operated at a proton energy of 
820 GeV and a positron energy of 27.5 GeV. 



A detailed description of the ZEUS detector can be found elsewhere ||11|| . A brief outline 
of the components most relevant for this analysis is given below. 

Charged particles are tracked by the central tracking detector (CTD) []T2j| , which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift-chamber layers, organised in 9 superlayers covering the polar-angleQ 
region 15° < 9 < 164°. The relative transverse-momentum resolution for full-length tracks 
is a{pt)/pt = 0.0058p t © 0.0065 © 0.0014/p t , with p t in GeV. 



1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton-beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as r\ = — ln(tan|), where the polar angle, 9, is measured with respect to the proton-beam 
direction. 
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The high-resolution uranium-scintillator calorimeter (CAL) |TB[ consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in FCAL and BCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter is called a cell. The CAL relative energy 
resolutions, as measured under test-beam conditions, are cr(E) /E — 0.18/v^E for electrons 
and a(E)/E = 0.35/v 7 ^ for hadrons (E in GeV). 

Low-Q 2 events (0.17 < Q 2 < 0.70 GeV 2 ) were tagged by requiring the identification of the 
scattered positron in the beam pipe calorimeter (BPC) fl4], I5|. The BPC was a tungsten- 



scintillator sampling calorimeter, located close to the beam pipe, 3 m downstream of the 
interaction point in the positron-beam direction. The relative energy resolution from 
test-beam results was cr(E) / E = 0.17/ \[E [E in GeV). Each scintillator layer consisted of 
8 mm-wide strips. Using the logarithmically weighted shower position, the impact position 
of the scattered positron could be measured with an accuracy of about 1 mm. For events 
with Q 2 > 3 GeV 2 , the impact point of the scattered positron was determined with the 
small-angle rear tracking detector (SRTD) |16j or the CAL. The SRTD is attached to 
the front face of the RCAL and consists of two planes of scintillator strips, 1 cm wide 
and 0.5 cm thick, arranged in orthogonal orientations and read out via optical fibres and 
photomultiplier tubes. It covers a region of about 68 x 68 cm 2 in X and Y, excluding a 
10 x 20 cm 2 hole at the centre for the beam pipe. 



The LPS (T^j detected positively charged particles scattered at small angles and carrying a 
substantial fraction, xl, of the incoming proton momentum; these particles remain in the 
beam pipe and their trajectory was measured by a system of silicon microstrip detectors 
that can be inserted very close (typically a few mm) to the proton beam. The detectors 
were grouped in six stations, SI to S6, placed along the beam line in the direction of 
the proton beam, between 23.8 m and 90.0 m from the interaction point. The track 
deflections induced by the magnets of the proton beam-line allow a momentum analysis 
of the scattered proton. For the present measurements, only stations S4, S5 and S6 
were used. The resolutions were better than 1% on the longitudinal momentum and 
5 MeV on the transverse momentum. The effective transverse-momentum resolution is, 
however, dominated by the intrinsic transverse-momentum spread of the proton beam at 
the interaction point, which is about 40 MeV in the horizontal plane and about 100 MeV 
in the vertical plane. 



4 Reconstruction of the kinematic variables 

Different methods have been used for the reconstruction of the kinematic variables Q 2 and 
W, depending on the Q 2 range of the measurement. At low Q 2 , 0.17 < Q 2 < 0.70 GeV 2 
(hereafter referred to as the "BPC region"), the energy, E' e , and angle, 8 e , of the scattered 
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positron measured in the BPC were used ("electron method") to determine the kinematic 
variables from 



Q 2 = 2E e E' e (l + cos6 e 



W 



4:E e E p 



2E P 



1 — COS6L 



where E p and E e represent the proton and positron beam energies, respectively. For 
Q 2 > 3 GeV 2 (the "DIS region"), Q 2 and W were reconstructed with the double angle 



method [18] using the energy depositions in the CAL. 



For the reconstruction of the mass of the diffractive system X, the energy deposits in the 
CAL and the track momenta measured in the CTD were clustered to form energy-flow 
objects (EFOs) [0, |l9j- The EFOs thus include the information from both neutral and 
charged particles in an optimal way. The mass, Mx, was then obtained from the EFOs 



via 



Ms 



\ 



5> - £ 



PXi 



£ 



PYi 



£ 



PZi 



where the subscript i denotes an individual EFO; the EFOs associated with the scattered 
positron are excluded from the sums. 

The momentum of those scattered protons detected in the LPS, p hPS , was measured, 
along with its component perpendicular (parallel) to the beam direction, p^ PS (p z PS )- 
From these quantities, the fractional momentum of the scattered proton, xl, and t were 
determined via 



x L = Vz/E pi 

+ G^ PS ) 2 



Two quantities, y and 5 = — Pz)i + E e i{l — cos^e), the sum of E — Pz over 

all final-state particles in the event, were used in the event selection. The former was 
reconstructed either using the electron method (and denoted by y e ) or from the EFOs 
using the Jacquet-Blondel estimator [E0J as 



VJB ~ 2E e ' 

where the sum is over all EFOs, excluding those assigned to the scattered positron. En- 
ergy and momentum conservation require 5 to be twice the positron beam energy for a 
completely measured final state and neglecting resolution effects. 
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5 Monte Carlo simulation 



Monte Carlo (MC) generators were used to determine the acceptance of the apparatus. 
The reaction ep — > eXp was simulated in the BPC region with the EPSOFT2.0 |7, 21],|22 



MC generator interfaced to HERACLES4.6 [23], which simulates initial- and final-state 



QED radiation. For the description of the diffractive dissociation of virtual photons, 
7*p — > Xp, EPSOFT uses the triple- Regge formalism [|IJ, in which the inclusive diffractive 
cross section can be expressed in terms of three trajectories. If all the trajectories are 
Pomerons (IP IP IP), the cross-section da jdM\ is approximately proportional to \jM\. If 
one of the trajectories is a Reggeon (IPIPIR), the cross-section da jdM\ falls as ~ \ jM\. 
EPSOFT also simulates exclusive vector-meson production, ep —>■ eVp, where V = p°, 
it) or <p, and non-diffractive ep interactions, ep — > eY. Production of J/ip mesons has 
negligible effects on the acceptance and was not considered. EPSOFT was also used 
to simulate the double-dissociative reaction, ep — > eXN, where the proton diffractively 
dissociates into the state N. 



The second generator, used for the DIS region, was RAPGAP2.06 [24|, where, for the 



diffractive structure function, a factorisable expression was assumed based on the model 



of Ingelman and Schlein | 2"5| . In particular, a superposition of non-interfering Pomeron 
and sub-leading trajectories was used ("fit B", as determined by the HI Collaboration 0) 
together with the "fit 3" Pomeron parton density functions 0. Again, initial- and final- 
state QED radiation were simulated using HERACLES. 

All generated events were passed through the standard ZEUS detector simulation, based 
on the GEANT3.13 program p6|, and the trigger-simulation package. 



6 Identification of the scattered positron 

For the BPC sample, the events were selected in the trigger by requiring the presence of a 
scattered positron in the BPC. A positron with energy greater than 7 GeV was required in 
the offline analysis [p^,[15|. The following cuts were applied to reduce the contamination 
from photoproduction events, radiative events, and beam-related background: 

• y JB > 0.05; 

• 35 < 5 < 65 GeV; 

• I-ZvtxI < 50 cm, where Zy^x is the Z coordinate of the reconstructed vertex. 

For the DIS sample, the events used for the analysis were selected in the trigger by 
requiring the presence of a scattered positron in the CAL. Offline, a positron in the RCAL 
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with energy greater than 10 GeV was required. A positron finder based on a neural- 
network was used |27| . The following cuts were applied to reduce the contamination from 
photoproduction events, radiative events, and beam-related background: 



• y JB > 0.03; 

• y e < 0.95; 

• 35 < 5 < 65 GeV; 

• -50 < Z yTX < 100 cm. 

7 The LPS method 

Diffractive events are characterised by a final state proton scattered at very small angle and 
with energy nearly equal to that of the incoming proton. In the LPS method, diffractive 
events are then defined as those having a proton detected in the LPS with ~ 1 ■ 
Figure |I](a) shows the measured xl spectrum, uncorrected for acceptance. The diffractive 
peak is clearly visible for values of xl close to unity. For the present analysis, xl > 0.97 



was required. Previous studies |28[ indicate that the double-dissociative contribution to 
such events is negligible. 

Two data samples, collected in 1995, were analysed with the LPS method. The BPC 
sample, corresponding to a luminosity of 0.90 ± 0.01 pb _1 , covers the range 0.17 < Q 2 < 
0.70 GeV 2 and 90 < W < 250 GeV. The DIS sample covers the region 3 < Q 2 < 80 GeV 2 
and 80 < W < 250 GeV, and corresponds to a luminosity of 3.30±0.03 pb _1 . The analysis 
was limited to the range 3 < Mx < 38 GeV for the BPC sample and 3 < Mx < 33 GeV 
for the DIS sample. 

The candidate proton was tracked along the proton beam line and was rejected if, at any 
point, the reconstructed minimum distance of approach to the beam pipe, A P i po , was less 
than 400 fim or if the distance to the edge of the sensitive region of any LPS station, 
Api an e, was smaller than 200 /xm. These cuts reduce the sensitivity of the acceptance 
to the uncertainty in the position of the beam-pipe apertures and of the detector edges. 
In addition, t was required to be in the region 0.073 < \t\ < 0.4 GeV 2 , where the LPS 
acceptance is well understood f[? |. Beam-halo background results from a scattered proton, 



with energy close to that of the beam, originating from an interaction of a beam proton 
with the residual gas or with the beam collimators. In this case, the proton measured in the 
LPS is uncorrelated with the activity in the central ZEUS detector. This background was 
suppressed by the requirement that the sum of the energy and the longitudinal component 
of the total momentum measured in the CAL, the BPC and the LPS be less than the 
kinematic limit of 2E p : (E + Pz)cal + (E + Pz)bpc + 2p^ PS < 1665 GeV. This cut takes 
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into account the resolution of the measurement of p^ PS . A residual beam-halo background 
of approximately 8%, remaining after the cut, was subtracted statistically. 

In the BPC analysis, the measured number of diffractive events was corrected bin-by-bin. 
From this acceptance-corrected number of events, the cross section for the diffractive dis- 
sociation of virtual photons, 7*p — > Xp, was determined taking into account the integrated 
luminosity, bin widths, and bin-centring corrections. 

In the DIS analysis, the cross section for the diffractive dissociation of virtual photons 
at a given point within a bin was obtained from the ratio of the measured number of 
diffractive events to the number of events in that bin predicted from the MC simulation, 
multiplied by the j*p — > Xp cross section calculated at that point by the Monte Carlo 
generator. Both the acceptance and the bin-centring corrections were thus taken from 
the MC simulation. 

In both the BPC and DIS analyses, the diffractive cross-section daj^ / dMx was directly 
measured only in the region 0.073 < \t\ < 0.4 GeV 2 and extrapolated to the full t range 
using the t dependence assumed in the Monte Carlo generator. In the region covered by 
the present measurements, this is roughly equivalent to carrying out an integration over 
t assuming an exponential dependence on t, e _b '*', with b ~ 7.5 GeV -2 . 



8 The M\ method 

Diffractive photon dissociation, 7*p — > Xp, is characterised by the exchange of a colour- 
less object, the Pomeron, between the virtual photon and the proton. This suppresses 
QCD radiation, and hence the production of hadrons, in the rapidity region between the 
hadronic system X and the scattered proton, yielding a forward rapidity gap, a charac- 
teristic feature of diffractive interactions. This feature is reflected in the dependence of 
the cross section on Mx, daj^ / dMx oc 1/M x ° l ' p ^''~ 1 , i.e. approximately flat as a function 
of lnM x . In contrast, for non-diffractive events, large rapidity gaps are exponentially 
suppressed by QCD radiation, which populates the region between the struck quark and 
the coloured proton remnant. In this case, under the assumption of uniform, random 
and uncorrelated particle emission in rapidity, the \nM x distribution falls exponentially 
towards low Mx values. The different properties of the In M x distribution for diffractive 
and non-diffractive events are exploited in the M\ method 0j. 

The M\ method was used to analyse BPC data taken in 1996-97, corresponding to an 
integrated luminosity of 6.2 ±0.1 pb _1 . The kinematic range used was 0.22 < Q 2 < 
0.70 GeV 2 , 90 < W < 220 GeV and 3.0 < M x < 12.2 GeV. 
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8.1 Selection of the diffractive signal 



Figure 0(b) shows a representative distribution of In for data-in the bin 0.220 < Q 2 < 
0.324 GeV 2 and 150 < W < 180 GeV, compared to the distribution of the simulated events 
generated with EPSOFT. Also shown are the four individual contributions generated with 
EPSOFT for non-diffractive events, for the JPJPJP and IPIPIR contributions (shown 
combined in the figure), which lead to the diffractive dissociation of the photon, and 
for vector-meson production. Diffractive events dominate the region of low lnM|, while 
non-diffractive events exhibit a large peak at high hiM x and a steep exponential fall-off 
towards lower \nM x values. The relative weights of the four subprocesses are obtained 
from fits to the In M\ distribution of the data. The resulting sum of the MC events (open 
histogram) from the various subprocesses provides a reasonably good description of the 
data in the region of interest, lnM| < 8.5. 

In the region lnM|- ~ 4, the diffractive contribution to the events in Fig. |I](b) depends 
only weakly on In M\ . The expression 

i^wr D+CeMBham (4) 

was therefore fitted to the data and the parameters D, C and B were determined for each 
(Q 2 ,W) region. The exponential term (with B = 1.44 ±0.02), ascribed to non-diffractive 
events, was subtracted statistically to obtain the diffractive contribution. The parameter 
D was thus not used directly to determine the diffractive cross section. The exponential 
term resulting from the fit to the data of Fig. |l](b) is shown in Fig. |TJ(c) . 

The cross-section measurement was restricted to the range 2.2 < lnM^ < 5.0, corre- 
sponding to 3.0 < Mx < 12.2 GeV. The lower limit on Mx suppresses the contribution 
from diffractive vector-meson production, while the upper bound was chosen such that 
the non-diffractive contribution to the higher-Mx bins was always less than 50%. 



8.2 Proton-dissociative contribution 

The diffractive sample of ep — > eXp events selected with the M\ method as discussed 
in Section |H] contains a contribution from the double-dissociative reaction ep — *> eXN. 
The system N escapes undetected through the forward beam pipe, unless the proton 
dissociates into a state of sufficiently high mass, in which case some of the particles from 
the system N have transverse momenta large enough that they are detected in the FCAL 
region around the forward beam pipe. The contribution of the double-dissociative reaction 
ep — > eXN was simulated and studied with EPSOFT. 

Energy deposits in the FCAL, arising from the proton-dissociative remnant, give rise to 
a measured value of Mx considerably higher than the true value. In such events, there 



The data shown in Figs. [y(b,c) result from the cuts discussed in Section 8.2 
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is a gap in rapidity between the FCAL deposits from the proton remnant at high r\ and 
the system X at lower rj, and the invariant mass of the low-77 system is small with respect 
to the measured (apparent) Mx- Exploiting these features, events were rejected if they 
fulfilled all of the following three conditions: 

• the maximum 77 (t/max) of the EFOs was greater than 2.5; 

• the maximum rapidity gap between adjacent EFOs was greater than 3.5; 

• the mass reconstructed from EFOs with 77 < 2.5 was less than 0.6Mx- 



These cuts rejected approximately 10% of the data sample. The simulations using the 
EPSOFT MC program indicate that about 45% of these rejected events are from the 
double-dissociative reaction, / ~y*p — ► XN, and the events that survive the cuts consist 
of photon-dissociative events as well as events from the double-dissociative reaction with 
M N < 6 GeV. 

The measured number of events in each (Q 2 , W, Mx) bin was corrected for acceptance to 
determine the number of produced events by means of the Singular Value Decomposition 



(SVD) method | 29fl , which allows the evaluation of error correlations between adjacent 
bins. The number of events thus obtained was divided by the luminosity and the bin- 
widths to evaluate the average ep — > eXp three-fold differential cross section for each 
(Q 2 , W, Mx) bin. From this, the cross section was obtained using Eq. ([[]), integrated 
over t and evaluated at the logarithmic centre of the bin. The residual double-dissociative 
contribution in the data leads to an overestimation of the cross section for the diffractive 
dissociation of virtual photons, 7*p — ► Xp; it was evaluated using the LPS data and 



subtracted as discussed in Section 10. 



9 Systematic uncertainties 



The main sources of systematic uncertainty can be classified into three groups: the 
positron measurement in the BPC or the CAL; the measurement of the hadronic final 
state in the CAL; and the measurement of the scattered proton in the LPS: 



• measurement of the scattered positron: 

— for the BPC samples, the effects of the uncertainty in the absolute BPC energy 
scale [|TjJ (±0.5%), the positron-selection criteria and the alignment of the BPC 
result in an uncertainty in the cross section that is typically ±7% and always 
smaller than ±20%; 



10 



for the DIS sample, the 10 GeV cut on the scattered-positron energy was changed 
by ±2 GeV [p8"| |. The parameters of the neural-network positron finder were mod- 
ified. To check the acceptance at low Q 2 , which is determined by the positron 
position in the SRTD, the fiducial region around the impact point of the positron 
was changed. The resulting systematic uncertainty is typically ±7% and always 
smaller than ±25%. 



• for the determination of the uncertainties related to the hadronic final state for the 
M\ method, the effect of the uncertainty in the CAL energy scale (±2%) was studied 
and the parameters in the algorithm that forms EFOs were varied. The non-diffractive 
slope (B in Eq. (f|)) was varied between 1.42 and 1.46. It was checked that a differ- 
ent choice for the functional form of the diffractive contribution in Eq. (|4|) does not 
significantly affect the final number of diffractive events. The resulting uncertainty is 
typically ±8% and always smaller than ±12%; 

• the systematic uncertainties in the measurement of the scattered proton in the LPS 
were estimated as follows: 

— to estimate the sensitivity of the LPS acceptance to the uncertainties in the po- 
sitions of the beam-line apertures, the lower limits on the distance of closest ap- 
proach to any of the beam-line elements and to the edge of the sensitive region of 
each detector were raised from 400 jum to 1000 jum and from 200 /mi to 300 /im, 
respectively; 

— the xl window was varied by ±0.01; 

— the uncertainty in the subtraction of the beam-halo events was estimated by re- 
moving the E + pz cut. 

The resulting systematic uncertainty arising from the LPS measurement is typically 
±10% and always smaller than ±25%. 



In addition, the Mx, W and t dependences in EPSOFT and RAPGAP were varied within 
the limits allowed by the data, yielding changes in the cross section negligible with respect 
to all other uncertainties. The relative fraction of vector meson production in EPSOFT 
was varied by up to ±10%, again with negligible effects on the results. 

All the above contributions were summed in quadrature to give the final systematic un- 
certainties. The normalisation uncertainty due to the luminosity determination is ±1% 
for the 1995 data and ±1.5% for the 1996-97 data and was not included in the sum. 

The value of R D was assumed to be zero throughout the analysis. Given the absence 
of experimental information on R D , no attempt was made to quantify the systematic 
uncertainty entailed by this assumption. 
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10 Cross-section measurements and comparison of 
the M\ and LPS methods 

The values of cPa^ /dMxdt and daJ^/dMx extracted with the LPS method are given 
in Tables [I] and 0, respectively. The results obtained with the M\ method are presented 
in Table |3|. All results are corrected to the Born level. 

As discussed in Section |872|, the sample from the M x method contains a double-dissociative 
contribution. Since the sample selected with the LPS method has a negligible proton- 
dissociative background ||28|| , the contamination in the BPC sample was estimated by 



directly comparing the results from the two analysis methods. 

To achieve this, the ratio, Rm x , of the average cross section measured with the M\ method 
and that measured with the LPS method was determined in a single kinematic region 
corresponding to the bins given in Table |3|. The value obtained, Rm x = 1.85±0.38 (stat.), 
is attributed to a substantial contribution from the double-dissociative reaction to the 
cross section measured with the M x method. In terms of the ratio of the number of double- 
dissociative events to the total number of events in the sample, i?diss = (1 ~ 1/Rm x ), the 
estimated value of Rm x corresponds to -Rdiss = (46 ± 11)%. This is consistent with a 
previous estimate at higher Q 2 of (31 ± 15)% |7|]. 

The results obtained with the M x method presented in this paper were corrected for the 
residual double-dissociative background using the measured value of Rm x • The correction 
was assumed to be independent of W and Q 2 , in agreement with the hypothesis of vertex 
factorisation [p0|1 . The values of daj^ / dMx extracted with the M\ method for each (Q 2 , 



W, Mx) bin are given in Table The subtraction of the double-dissociative background 
entails a ±21% uncertainty in the normalisation of the cross sections obtained with the 



M x method. 



11 Results and discussion on the W dependence of 
the diffractive and total cross sections 

The energy dependence of the photon-dissociative cross sections can be successfully de- 
scribed by a power of W, both for photoproduction |3],f|] and for DIS Jj,j7|] - at least in 
the region of small xjp values where the exchange of sub leading Regge trajectories can 
be neglected. Although the experimental uncertainties are large, the value of this power 
is different for the two regimes. This is analogous to the behaviour observed for the W 
dependence of the virtual photon-proton total cross section, er^ot" |3],[TT]: the slow rise 
of 

a tot a t high W observed in photoproduction becomes faster at high Q 2 ; the transition 
takes place for Q 2 ~ 1 GeV 2 . In this section, the W dependence of the photon-dissociative 
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cross section, daJ^/dMx, is studied in this transition region and is compared to the W 
dependence of o"7 ot p by considering the ratio of daj^ / dM x to cr^ ot p . 



11.1 The W dependence of the diffractive cross section 

Figure |2] shows the values of the diffractive cross sections extracted with the M\ method 
in the BPC region as a function of W for three Q 2 and two Mx bins. The form 

^JfL = Ai ■ W a ^ (5) 

was fitted to these data, where a^m is a global parameter and the normalisation parameters 
Ai were left free to vary for each (Q 2 ,Mx) bin. The results of the fit, taking into account 
the correlations between adjacent bins, are shown in Fig. |2|; they give a good description 
of the data (x 2 jndj = 0.51, calculated using the statistical uncertainties only). The fitted 
value of the power of W is 

a diS = 0.510 ± 0.043(stat.)+[];;°2(syst.). 

Expressing the W dependence of the cross section in terms of an effective Pomeron inter- 
cept |U|, ajp, as 

^ oc (ty 2 ) 2a -- 2 , 

dM x 

the fitted value of a diS corresponds to 

civ = 1.128 ± 0.011(stat.)+J;j3o(syst.). 
This value of ajp can, in turn, be related to the Pomeron intercept, ajp(0), via 

qjp = ajp(0) — a'jp ■ 

where |t| is the mean value of \t\. The value of ap(0), obtained assuming a' F = 0.25 GeV -2 
and using |t| = 0.13 GeV 2 P5[|3*2]| , is ajp(0) = 1.161 ± 0.011(stat.)^Q '^(syst.); it is shown 



in Fig. together with the values determined from photoproduction and from higher- 
Q 2 measurements [§, ^, |, 0| . The quoted systematic uncertainty does not include the 
uncertainty on a'jp, which was also not included in the other results presented in Fig. [| 
The value of op(0) from the present measurement at low Q 2 does not differ significantly 
from the values at higher Q 2 . 

Equation (H) was also fitted to the data allowing different values for the parameter a^s 
in the three Q 2 bins of the measurement; the three resulting values of a^m are compatible 
with the global value. The data used in the fit have values of xjp typically much smaller 
than 0.01, with the exception of the bin with lowest W and highest Mx values, which 
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receives contributions from xjp values up to xjp = 0.018. It was assumed that Pomeron ex- 
change dominates in this region, and no attempt was made to include secondary Reggeon 
exchange in the fit. Finally, it should be noted that a possible W dependence of the 
double-dissociative fraction would affect the extracted value of ajp. 

Figure ^jalso shows ajp(0) as obtained from the ALLM97 parameterisation [ 33 1 of the j*p 



total cross section, which gives a good representation of the inclusive F 2 data for the entire 
Q 2 range. The value of ajp(0) from ALLM97 is consistent with the present determination 
from the diffractive data in the BPC region, whereas in the DIS region it is higher than 
the HI and ZEUS diffractive measurements ||[7|]. 

The LPS cross sections are presented in Fig. f|; they are in agreement with the previous 
ZEUS measurements at large Q 2 and with the present BPC data obtained with the M x 
method. The previous ZEUS data obtained by the M\ method have also been corrected 
for the residual double-dissociative background using the value of Rm x given in Section [K| 
to make a direct comparison with the earlier data, the BPC cross sections from Table |3| 
have been interpolated to Mx = 5 and 11 GeV using bin-centring corrections based on 
EPSOFT. The solid lines in Fig. f| correspond to the fit of Eq. (^) to the BPC data alone, 
which also provides a good description of the DIS data (dashed lines). Figures § and f| 
thus show that the W dependence of the inclusive diffractive cross section exhibits no 
significant changes from the BPC to the DIS region. 



11.2 Comparison of the W dependence of the diffractive and the 
total cross sections 

The W dependences of the diffractive and total cross sections were directly compared by 
studying their ratio 

da^/dM x 

This ratio is plotted as a function of W in Fig. |j| where the values of the diffractive cross 
sections shown in Fig. |] were divided by the corresponding values of the j*p total cross 
section, [15, 53]. The lines denote the fit shown in Fig. |] divided by the corresponding 



values of the ALLM97 parameterisation for a"7 ot p . The lines give a good representation of 
all the data. While there is a clear increase in r as a function of W for Q 2 < 1 GeV 2 , for 
higher Q 2 the distribution is flat in W. 

The form r = iVj • W p was fitted to the BPC data measured with the M\ method; 
here, p is a global parameter and the normalisation parameters, Ni, were left free to 
vary for each (Q 2 , Mx) bin. The fit gives a good description (not shown) of the data 
with p = 0.24 ± 0.07, where the uncertainty is derived from the fit, consistent with the 
expectation [|7| from Regge theory that 
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(d^/dMx) (W 2 ) 2 ^ 2 = w ^ apW -i } „ w o. W 

{W 2 )a]p{0) -X 



This result suggests a different behaviour from that found in the DIS region, where the 
value p = 0.00 ± 0.03 J7J indicates that the diffractive and inclusive cross sections have 
the same W dependence, contrary to the expectations of Regge theory. 

In summary, in the BPC region the W dependence of the diffractive cross section is 
compatible with that expected from Regge phenomenology. The ratio between diffractive 
and total cross sections grows with W at a rate consistent with Regge theory. This is 
in contrast to the DIS region, where the expectations of Regge theory for the ratio of 
diffractive and total cross sections are not fulfilled, since the ratio is flat as a function 
of W. This difference in the W dependence of the ratio is reflected in the fact that the 
values of azp(0) extracted from the diffractive cross section and from are similar in 
the BPC region, but not in the DIS region. 



12 Results and discussion on the Q 2 dependence of 
diffractive and total cross sections 

The Q 2 dependence of cr^f has been observed to change around Q 2 ~ 1 GeV 2 com- 
pared to the approximate l/Q 2 scaling behaviour found at high Q 2 , data at Q 2 ^ 1 GeV 2 
exhibit a weaker Q 2 dependence, with being nearly independent of Q 2 at the lowest 
Q 2 values measured. This is consistent with the expectation from the conservation of the 
electromagnetic current that approaches a constant or, equivalently, that vanishes 
like Q 2 as Q 2 -> 0. 

In this section, the Q 2 dependence of the diffractive cross section, daJ-^/dMx, is studied 
and the question is addressed of whether and where a transition similar to that observed 
for &2ot occurs for the diffractive dissociation of virtual photons. 



12.1 The Q 2 dependence of the diffractive cross section 

Figure ^| shows the diffractive cross sections, daj^ /dM x , as a function of Q 2 in bins of W 
and Mx- The present measurements are plotted together with previous ZEUS results 
obtained with the M\ method in the DIS region, and HI results | 35[| , obtained with the 
rapidity-gap method in photoproduction for Mjv < 1.6 GeV and \t\ < 1 GeV 2 ; bin-centring 
corrections based on EPSOFT, analogous to those described in Section |Il.l| , were applied, 
where necessary, to both ZEUS and HI results. No further corrections were applied to the 
HI data; notably, no attempt was made to correct for the double-dissociative background. 
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In Fig. H a change in the Q 2 dependence of da^g/dMx as Q 2 increases is apparent and is 
similar to that observed in the o"7 ot p data: at low Q 2 , the data do not exhibit a strong Q 2 
dependence, while at larger Q 2 , the cross section falls rapidly for increasing Q 2 . Figure [7| 
shows function of Q 2 for fixed W and Mx', while at large Q 2 the data 

do not exhibit a strong Q 2 dependence, XpF®^ falls by a factor of about ten between 
Q 2 « 8 GeV 2 and Q 2 « 0.2 GeV 2 . 

12.2 Discussion 

The diffractive dissociation of the virtual photon can be described by perturbative QCD 
(pQCD) since the photon's virtuality, Q 2 , provides a hard scale. In particular, in the 
proton rest frame, the reaction can be viewed as a sequence of three successive processes || 
|36|| : the photon fluctuates into a qq (or qqg) state, the qq dipole scatters off the proton 
target and, finally, the scattered qq pair produces the final state. At high centre-of-mass 
energies of the 7p system, these processes are widely separated in time. The qq fluctuation 
is described in terms of the photon wave-function derived from QCD. The interaction of 
the qq dipole with the proton is mediated, in lowest order, by the exchange of two gluons 
in a colour-singlet state. 

The present results have been compared to the model of Bartels et al. (BEKW) |J7[] . In 
this model, neglecting the contribution of longitudinally polarised photons, the dominant 
(leading-twist) contributions to the diffractive structure function in the kinematic domain 
of the present measurement come from the fluctuations of the photon into either a qq pair 
(Fqq) or a state {.Fqq g )- The (3 spectra of these two components are determined 
by rather general properties of the photon wave-function: F gg behaves like (3 (1 — (3) 
and F ggg like (1 — (3) 1 , where 7 = 3.9 0, |37]]. At large (3, qq production dominates 
over qqg production, while, at small (3, qqg production becomes dominant. F gg has no 
Q 2 dependence; F^ gg is of order as and has a logarithmic Q 2 dependence of the type 
log (1 + Q 2 /Ql), where the scale parameter Ql is taken to be 1 GeV 2 . The model does 
not fix the xjp dependence of F^ and F% , but assumes for both a power-like behaviour 
x ^diff(Q ) ^ w ]j ere exponent is determined from fits to the data. 

A comparison of the BEKW parameterisation with the present data is shown in Figs. |6| 
and 0. The values of the parameters, including the normalisation of the F gg and 
components, were taken from a fit to the previous ZEUS results [[7|, with the exception 
of the xjp exponent, for which a constant value corresponding to determined from 

Eq. (|), was used. The DIS data at high Q 2 constrain the parameterisation of the (3 
dependence of F^ (dashed lines) at low Mx and of F^ g (dotted lines) at high Mx- The 
logarithmic Q 2 dependence of F^ gg is probed only in the highest-Q 2 region and is less well 
constrained. 

The Q 2 dependence of F ggg becomes crucial in the transition to low Q 2 . In fact, as Q 2 
decreases from the DIS into the BPC region, for a given value of Mx, (3 also decreases: 
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the BPC data in Figs. |^ and [7| correspond to values of (3 that are typically two orders of 
magnitude smaller than those in the DIS data and thus, in the BPC region, the contri- 
bution from the fluctuation of the photon into a qqg system becomes dominant. While 
extrapolating the BEKW parameterisation to low Q 2 may not be formally justified, it is 
interesting to note that, in F^ , conservation of the electromagnetic current is assured by 
the fact that log (1 + Q 2 /Qq) vanishes like Q 2 ' jQ\ as Q 2 — > 0. The transition from the 
linear behaviour at low Q 2 to the logarithmic behaviour at higher Q 2 is controlled by the 
scale parameter Ql; the choice Ql = 1 GeV 2 successfully describes the BPC data. 



12.3 Comparison of the Q 2 dependence of the diffractive and 
the total cross sections 

Figure |] shows the ratio r = (daJJ^ / dMx) / crtot as a function of Q 2 for different W and 
Mx bins. At low Q 2 , the Q 2 dependence of the diffractive cross section is similar to that 
of cr^ ot p . In the DIS regime, daj^ / dM x decreases with Q 2 more rapidly than a^. This is 
more evident for small values of Mx- In addition, the ratio r appears to increase between 
the BPC and the DIS regions. 

Also shown in the figure are the results of the BEKW parameterisation of the diffractive 
cross-section daJ^/dMx, shown in Fig. ||, divided by the values of cr^otf gi yen by the 



ALLM97 parameterisation ||33|| . There is reasonable agreement between these parame- 
terisations and the data, indicating that the data may be qualitatively described by an 
appropriate choice of the relative fractions of the qq and qqg contributions. 



13 Summary 

The diffractive dissociation of virtual photons, 7*p — > Xp, has been studied at HERA 
at low Q 2 (0.17 < Q 2 < 0.70 GeV 2 ) and in deep inelastic scattering (DIS) (3 < Q 2 < 
80 GeV 2 ). The diffractive signal has been selected either by requiring the detection of a 
final-state proton with at least 97% of the incoming proton-beam energy, or by exploiting 
the different properties of the Mx distributions for diffractive and non-diffractive events. 

The W dependence of the low-Q 2 cross-section data obtained with the M x method (3 < 
Mx < 12.2 GeV) has been found to be compatible with a single power of W, which 
corresponds to a Pomeron intercept of Ozp(0) = 1.161 ± 0.011(stat.)^QQ3Q(syst.). This is 
consistent with that previously observed in the DIS regime. Thus, the significant change 
in the W dependence exhibited by the 7*p total cross section in the transition from low 
Q 2 to DIS is not visible in the diffractive cross section. To elucidate this difference, the 
W dependence of the ratio, r, of the diffractive cross section to the 7*p total cross section 
was studied at low Q 2 and was found to rise with W, r oc jy°- 24±0 - 07 ; j n agreement with 
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the expectation from Regge theory. This is in contrast to the observation at higher Q 2 
that this ratio is independent of W. 

The Q 2 dependence of the diffractive cross section changes as Q 2 increases up to the DIS 
regime: while at low Q 2 the data do not exhibit a strong Q 2 dependence, at larger Q 2 the 
cross section falls rapidly for increasing Q 2 . This change of behaviour occurs for values 
of Q 2 around 1 GeV 2 and is analogous to that observed in the total 7*p cross section. 
The ratio of the diffractive cross section to the 7*p total cross section was studied as a 
function of Q 2 . At low Q 2 , the ratio r shows little dependence on Q 2 , indicating that 
the Q 2 dependence of the diffractive cross section is similar to that of a^. The ratio 
increases between the BPC and the DIS regions. In the DIS regime for low M x , the ratio 
decreases with increasing Q 2 , indicating that the diffractive cross section decreases with 
Q 2 more rapidly than the 7*p total cross section. 

The main features of the data, reproduced by a parameterisation based on the BEKW 
model, indicate that the framework in which the incoming virtual photon fluctuates into 
a quark-antiquark pair is, in general, adequate to describe diffractive processes in ep 
collisions from the BPC to the DIS region. At the same time, the data suggest the 
increasing importance of the contribution from qqg states at low Q 2 . It is interesting that 
the ratio of the diffractive cross section to the total cross section shows a change from a 
w o.24±om dependence for Q 2 < 0.7 GeV 2 to W om±om for Q 2 > 3 GeV 2 . This complex 
behaviour of diffraction as a function of both Q 2 and W reveals a rich testing ground for 
future theoretical models. 
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Table 1: The values of (Pa 1 ^ /dMxdt measured with the LPS method in the 
range 0.073 < \t\ < 0.40 GeV 2 with the bin ranges indicated. The data are at 
(t) = 0.17 GeV 2 . The first and second error values represent the statistical and 
systematic uncertainties, respectively. The normalisation uncertainty related to the 
luminosity measurement is not included in the systematic uncertainty. 
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Table 2: The values of daJ^/dMx measured with the LPS method with the 
bin ranges indicated. The first and second error values represent the statistical and 
systematic uncertainties, respectively. The normalisation uncertainty related to the 
luminosity measurement is not included in the systematic uncertainty. 
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Table 3: The diffractive cross-sections daj^ /dMx measured with the M x method 
with the bin ranges indicated. The first and second error values represent the statis- 
tical and systematic uncertainties, respectively. The ±21% systematic uncertainty 
due to the double-dissociation correction is not included in the systematic uncer- 
tainty, nor is the normalisation uncertainty related to the luminosity measurement. 
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Figure 1: (a) The xl spectrum as measured in BPC events with the LPS; 
(b) and (c) the Yd.M x distribution (Mx in GeV) of the BPC data in the region 
0.220 < Q 2 < 0.324 GeV 2 and 150 < W < 180 GeV. In (a), the position of the 
arrow indicates the value xl = 0.97 used in the selection. In (b), the data are 
compared to the mixture of four kinds of EPSOFT MC events described in the text: 
region A corresponds to non-diffractive events, B+C to the sum of the IP IP IP and 
IPIPIR contributions and D to the vector-meson contribution. In (c), the straight 
line shows the exponential slope, resulting from the fit described in the text, for 
non-diffractive events. 
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Figure 2: Diffractive cross sections for (a) Q 2 = 0.27 GeV 2 , (b) Q 2 = 0.39 GeV 2 , 
and (c) Q 2 = 0.58 GeV 2 for two different M x ranges as a function of W . The 
inner bars indicate the size of the statistical uncertainties; the outer bars show 
the size of the statistical and systematic uncertainties added in quadrature. The 
points were corrected for the double- dissociative background; the associated ±21% 
normalisation uncertainty is not included, but is shown separately as a shaded error 
band. The normalisation uncertainty associated with the luminosity measurement 
is not shown. The lines show the results of the fit described in the text. 
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Figure 3: Results for a jp(0) m different Q 2 regions. The value ofacjp(0) obtained 
from this analysis is shown as the solid circle. The open symbols show the results 
from the photoproduction /j^,^/ and DIS diffractive analyses The inner bars 

indicate the size of the statistical uncertainties; the outer bars show the size of the 
statistical and systematic uncertainties added in quadrature. The line is from the 
ALLM97 parameterisation Jffi] of the 7*p total cross-section data. 
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Figure 4: Diffractive cross sections for different Q 2 and Mx values as a function 
of W . The results obtained with the LPS method are shown as stars. The inner 
bars indicate the size of the statistical uncertainties; the outer bars show the size 
of the statistical and systematic uncertainties added in quadrature. The low-Q 2 
points obtained with the M\ method (full squares) were corrected for the double- 
dissociative background; the corresponding ±21% normalisation uncertainty is not 
included, but is shown separately as a shaded band. The normalisation uncertainty 
associated with the luminosity measurement is not shown. The open squares at 
high Q 2 are from a previous ZEUS publication and have been corrected for the 
double-dissociative background using the same estimate as for the low-Q 2 points, 
as discussed in the text. The solid lines are the results of the fit to the BPC data 
described in the text, which also gives a good representation of the higher-Q 2 data 
(dashed lines). 
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Figure 5: The ratio of the diffractive cross section to the^y*p total cross section for 
different Q 2 and Mx values as a function ofW. Other details are as in the caption 
to Fig. [|. The lines denote the fit shown in Fig. [| divided by the corresponding 
values of from the ALLM97 parameterisation. 
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Figure 6: ITie values of daj^ / dMx for different W and Mx values as a function 
of Q 2 . Other details are as given in the caption to Fig. [|. The solid lines are the 
results of the BEKW parameterisation described in the text; the dotted (dashed) 
lines are the results of the same parameterisation for the qqg (qq) contribution 
alone. Note the break in the Q 2 scale below ~ 10~ 2 GeV 2 . 
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Figure 7: 27ie values of xpF 2 for different W and Mx values as a function 
of Q 2 . Other details are as given in the caption to Fig. 0. The solid lines are the 
results of the BEKW parameterisation described in the text; the dotted (dashed) 
lines are the results of the same parameterisation for the qqg (qq) contribution 
alone. 
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Figure 8: The ratio of the diffractive cross section to the 7*p total cross section 
for different W and Mx values as a function of Q 2 . Other details are as given in 
the caption to Fig. 0. The solid lines are the results of the BEKW parameterisation 
described in the text, divided by the corresponding cr^ ot p values from the ALLM97 
parameterisation; the dotted (dashed) lines are the results of the same parameter- 
isations for the qqg (qq) contribution alone. Note the break in the Q 2 scale below 
~ 10" 2 GeV 2 . 
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